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Abstract 

A general formalism is worked out for the description of one-dimensional 
scattering in non-hermitian quantum mechanics and constraints on trans- 
mission and reflection coefficients are derived in the cases of V, T or PT in- 
variance of the Hamiltonian. Applications to some solvable "PT-symmetric 
potentials are shown in detail. 

Our main original results concern the association of reflectionless poten- 
tials with asymptotic exact VT symmetry and the peculiarities of separable 
kernels of non-local potentials in connection with Hermiticity, T invariance 
and VT invariance. 
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1 Introduction 



Since the seminal paper by Bender and Boettcher pQ, research on PT-symmetric 
quantum mechanics has been mainly focused on bound states, either with real 
energies, or in complex conjugate energy pairs, while relatively few authors have 
studied scattering states of Hamiltonians with both discrete and continuous spectra[21 

This gap has been recently bridged, at least in part, by a review paper on 
complex absorbing potentials [Hj, covering the more extended topics of one- 
dimensional scattering in non-hermitian quantum mechanics, with some general 
formulae valid for complex PT-symmetric potentials. Nonetheless, we think it 
worthwhile to go into further details in the latter case, while paying the due credit 
to the authors of Ref . jE] , and apply our formalism to some examples of solvable 
potentials for the sake of clarity. 

In any case, we share the philosophy of ref. jE], i. e. to consider PT-symmetric 
potentials as complex potentials defined on the real axis. Therefore, when we 
define complex coordinate shifts, we mean this as a method of generating complex 
potentials depending on a real coordinate. 

The main purpose of the present paper is thus to give a general description of 
one-dimensional scattering in non-hermitian quantum mechanics with the Hamil- 
tonian, H, invariant either under parity, P, or time reversal, T, or their product, 
VT, in the case H is not separately invariant under P and T. We do not intend, 
however, to provide here physical motivations for VT symmetry, nor to elaborate 
on the significance of violation of hermiticity and unitarity in quantum mechan- 
ics. The reader interested in these topics can consult a wide list of references, in 
particular the recent works |H1 E] and references therein. 

In the present work, we limit ourselves to a framework in which the scat- 
tering potentials vanish asymptotically: in particular, we exclude potentials di- 
verging at infinity, like those considered in Ref.|lUj. Therefore, our asymptotic 
wave functions are linear superpositions of plane waves and the present approach 
closely resembles standard scattering theory described in text-books on quan- 
tum mechanics, such as Ref. jTT], in the limit in which Hermiticity and P— and 
T— symmetries hold. 

What we try to achieve is to assemble a comprehensive and self-contained for- 
malism for discussing scattering problems in one-dimensional quantum mechan- 
ics. With this formalism available, we provide suitable examples, which should 
make the reader capable of elaborating his/her judgement on the relevance of the 
subject of VT symmetry, including some contributions to a better understanding 
of "exact" VT symmetry. Due to the relation of the Schrodinger equation to the 
classical Helmholtz equation, our formalism may be accommodated to deal with 
optics [SJ E2] with refraction index characterized by "handedness" [T3] . 

In addition, we try to clarify the interplay between T invariance and her- 
miticity, displaying a complex non-local solvable potential in which hermiticity 
does not force T invariance, but can be compatible with VT invariance. To our 
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knowledge, a VT— symmetric non-local potential is introduced and worked out 
here for the first time. 

We hope that our work is sufficiently self-contained, such as not to require 
any particular specific background of the reader. The paper definitely has topical 
review aspects, though it does not pretend to give a complete list of references. 
There are, however, relevant original results, in particular for non-local separable 
potentials. 

The paper is organized as follows: section El describes the basic formalism 
of one-dimensional scattering in non-hermitian quantum mechanics, section |3] 
introduces symmetries under which H may be invariant, section 0] defines density 
currents and continuity equations, section |5] applies the formalism worked out 
in the previous three sections to some solvable potentials. Finally, section El is 
devoted to conclusions. 



2 L-R Representation 

We start from the general time- dependent Schrodinger equation 

d 2 f d 

~dx 2 ~^ ^' ^ + / K (yX ' V > ^ ^' ^ dy = % ~dt^ ^' l ' ' ^ 

written in units h = 2m = 1. For a monochromatic wave, of energy u, the time 
dependence of the wave function is 

4>(x,t) = ty(x)e- iujt . (2) 

In the present work, unless otherwise stated, we consider local potentials, for 
which the kernel K reduces to 

K(x,y)=6(x-y)V(y) (3) 

If eqs.(J2HHl) hold, eq. (JT} reduces to the time independent Schrodinger equation 
satisfied by ty(x) 

= (--jL + V(x)j * = k 2 ^/, (4) 

with k = y/ui (> 0) the wave number. In order to solve eq.(JI]), it is convenient 
to work in a two dimensional Hilbert space where the basis vectors are the kets 
\R > and \L > (and the corresponding bras < R\ and < L\). In configuration 
space,with the choice of the time dependent phase given by eq. (J2J), 

<x\R,k >~ e ikx (5) 

represents a plane wave travelling from left to right and 

< x\L,k >~ e~ ikx (6) 
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a wave travelling from right to left. In the following, the explicit k dependence of 
the basis vectors will be omitted, whenever not strictly necessary, for simplicity 
of notation. 

In coordinate space, the asymptotic states , i.e., outside the (finite) range of 
the potential, can be expressed at x — > — oo as 

I**-- oo> = A.\R>+B.\L> (7) 

and, at x — > +00, as 

\* x -++ 0O ) = A+\R>+B+\L>, (8) 

or, in terms of wave functions 

j / \ / A^e ikx + B_e~ ikx x — > -00 , . 

^ ~ \ A + e ifcx + x -> +00 ^ ' 

In the case of a finite-range local potential, Eq. (jlj admits a general solution 
written as a linear combination of two independent solutions, Fi (x) and F 2 (x), 
with non-zero Wronskian, whose asymptotic expressions are both of the form: 

lim F m (x) = a m ±e ihx + b m± e~ ikx , (m = l,2) (10) 

x— >±oo 

The a m ± and b m ± are simply related to the asymptotic amplitudes A± and B± 

A± = aai± + (3a 2 ± (11) 
B± = abi± + (3b 2 ± ■ 

If *&i(x) = aFi (x) + (3F 2 (x) , inserted into Eq. (J2J), gives rise to a wave moving 
from x = —00 to x = +00, the amplitude of the regressive wave vanishes at +00: 

a b 2+ 

and the transmission and reflection coefficients of the wave moving from left to 
right are immediately written as 

_ A+ aa 1+ + (3a 2+ 
~* R A_ aai_ + /3a 2 _ 

a 2+b\+ — di+b 2+ 
a 2 _6i + - a x _b 2+ 



Rl^r = — — 



B_ 
Z 



+ /36 2 - 
aai_ + /3a 2 _ 

a 2 _6i+ — ai-b 2+ 
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The asymptotic form of ^i(x) is thus, neglecting a global normalization factor 

~ e tkx + R L ^ R e~ ikx , x^-oo (12) 
~ T L ^ R e ikx , x->+oc. 

Similarly, if fy 2 (x) gives rise to a wave moving from x = +00 to x — —00, the 
amplitude of the progressive wave vanishes at —00: 

a a,2- 

and the transmission and reflection coefficients of the wave moving from right to 
left are: 

„ _B- + (3b 2 - 



B+ ocb 1+ + 



J2+ 



a 2 -bi_ — di-b 



2- 



a 2 -b 1+ - ai_b 2+ 
_ A + aa 1+ + (3a 2+ 

Rr^L = — — ; —~T 7 

B + ab 1+ + (3b 2+ 
0,1+0,2- ~~ 



a 2 -6i+ — ai-&2+ 
As a consequence, the asymptotic form of ^ 2 (x) is 

* 2 (a0 ~ T^ L e"^, (13) 
~ e" ife:r + R R ^ L e ikx , x -> +00 

When we compute the Wronskian of ^1 and \l/ 2 , defined as 

d d 
= * 2 (a:) - * 2 (a:)— 

we readily obtain W^(— 00) = —2ikT R ^ L and W(+oo) = —2ikT L ^ R . Therefore, 
a necessary condition for the Wronskian to be constant on the x axis is T L ^ R = 
T r ^l. It is easy to check that dW/dx = for any well-behaved local potential. 
Therefore, the equality of the two transmission coefficients is satisfied for any 
such potential. 

The scattering matrix, S, connects the outgoing states at t — > +00 to the 
ingoing ones at t — > — 00 

\* ut) = S\* in ) . (14) 

The S matrix elements are directly linked to the transmission and reflection 
coefficients. In this basis, \R > and \L > can be rewritten as 



and 



If we have an ingoing wave of \R) type, then the outgoing wave is 



A± \ „( A_\ ( T L ^ R \ „ ( 1 



D- I S \0 )^\ RlZr I ~ S \ 
If the ingoing wave is of \L) type, we obtain, with the same procedure 

Rr^l \ ~ f 



T, 



R-+L 



SI I). (15) 



As a consequence of the above equations, the S matrix elements are explicitly 
given by 

Srr Srl \ / Tl^r Rr-^l 



S=l""« ^ = VT" " . (16) 

Note that our Ti^j(Ri^j) corresponds to T l (R l ) of Ref. jB] and to tj(rj) of 
Ref . . Our definition of S matrix is the same as that of Refs. [BJ Ej , while 
in the S matrix defined in the book by Merzbacher ^Tj the rows are exchanged 
with respect to ours. The latter author introduces also the transfer matrix, M, 
yielding the asymptotic states at x — > — oo when applied to those at x — > +oo 
(see Ref.|llj. formula (6.24)). In our basis, M, like S, is a 2 x 2 matrix, whose 
elements are easily expressed in terms of right and left transmission and reflection 
coefficients. 

4 ^ M( i + )■ (17) 



B - J V B i 

If the incident wave is of \R) type, the effect of M is expressed by the equation 

1 \ =M ( Tl - r 

Rl^r J V 

while, if the incident wave is of \L) type, the following equation holds 

From the equations given above, we obtain 

M RR M RL \ = ( 1/T L ^ R -R R ^ L /T L ^ R \ 

M^r Mll J \ Rl^r/Tl^r T r ^l — Rr->lRl->r/Tl^,r J 

with det M = T R ^ L /T L ^R. 

In general, any matrix, O, in the R — L basis 
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can be written also as a linear combination of basic dyadic operators 

nij = \i><j\, J = R,L (20) 

in the form 

o = o RR n RR + o RL n RL + o LR n LR + o LL n LL (21) 

This notation will be used in the following sections in the discussion of the 
invariance of the Hamiltonian with respect to various transformations. 



3 V, T and VT Symmetries 

In the present section, we study the transformation properties of the Hamilto- 
nian with respect to V, T and VT reflections. With particular reference to VT 
transformations, it is useful to introduce a discussion of the behaviour of vectors 
and matrices in the R — L basis under coordinate shifts. 



3.1 Coordinate Shifts 

In connection with the coordinate shift x — * x + X , with X a real number, 

let us define a displacement operator, V(X ), through its action on the basis 
vectors \R) and \L) 

V(X ) \R) = e ikX ° \R) ; 
V(X )\L) = e- ikX °\L); 

and on their dual vectors 

(R\V-\X ) = e- ikXo (R\; 
(L\V-\X ) = e ikXo (L\. 



In matrix form 



ikX 

W=( e -lx )• (22) 



Since X is real, P _1 (X ) = V*(X ). 

The basic dyadic operators are thus transformed under D according to the 
obvious relations 

V(X )\R)(R\V-\X ) = \R)(R\; 

V(X )\L)(L\V'\X ) = \L){L\; 

V(X )\R)(L\V-\X ) = e 2 * kXo \R)(L\; 

V(X )\L)(R\V- 1 (X ) = e- 2ikXo \L)(R\. 
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and a generic one-body operator, conveniently written in the form of a 2 x 2 
matrix, is consequently transformed as follows 



o^y^H^o^ C 0rl 

(23) 

3.2 Parity 

We shall consider now the parity transformation 

X ► X, p x ► Px • 

Under parity, the kets \R > and \L > are transformed according to 

V\R >= \L >; 
V\L >= \R >, 

and the bras < R\ and < L\ are changed as follows 

< RIV- 1 =< L\; 

< L\V- X =< R\. ' 

Hence, the four basic operators fi^ transformed according to VflijV -1 yield 

(24) 



Qrr Qrl ^LL &LR. 



&LR &LL &RL &RR 

In this basis, the parity operator can be represented by the matrix V ( with 
V 2 = 1) [H], which, acting on the left, exchanges lines, while acting on the right 
exchanges columns 

V=(° l)=V-* (25) 



x 1 , 

The Hamiltonian can be expressed as 

H = H RR Q RR + H RL Vt RL + H LR Vt LR + H LL {l LL , (26) 
or, in matrix form 

it _ ( Hrr Hrl 
V H LR H LL 

The transformed Hamiltonian is then 

Hp = VHP- 1 

= H RR Vl LL + H RL Vt LR + H LR Vl RL + H LL Vl RR 
Hll H lr 
Hrl H rr 



Parity invariance for the Hamiltonian, H-p = H, therefore requires 

Hrr = H LL (27) 



Hrl = H LR . (28) 

In the interaction picture, the S matrix is known to be the following limit of 
the transition operator, T(t — t ) (see Ref.|llj. formulae (14.49) and (20.7) ) 

S = lim f{t - t ) = lim U (t) T(t - t^ 1 (t ) (29) 

t — s-OO ,£{)—>■ — OO t — s-00,t{)— >— oo 

with 

U Q (t) = e iHot , H =p 2 = -^ (30) 

and 

oo 

T(t - t ) = e - iH{t ~ to) = ^{-i) n H n {t - to) n /n\ (31) 

n=0 

It is now possible to investigate the transformation of S under parity V. In 
order to compute 

S v = VSV~ l (32) 

we need 

oo 

VT(t-to)V~ x = ^2{-i) n VH n V~\t-to) n /n\ 

n=0 

oo 

= ^(-^) n ^(t-to)7n! 

n=0 

oo 

n=0 

and PUo (t) V' 1 = U (t), so that 

S v = S. (33) 

Hence, the S matrix commutes with V. 

Therefore, if the Hamiltonian is invariant under parity, so is the S matrix, 
which yields, upon using eq. (fTTIj) 

Srr = Sll — ► T L _> R = T R ^ L (34) 

and 

Srl = Slr — > Rr^l = Rl->r (35) 

The last two equations are strictly analogous to those, eqs. ()27|) and (|2*5j). 
expressing the parity invariance of the Hamiltonian, and are already given in 
Ref. 0. 



3.3 Generalized Parity 

We shall consider the following generalized parity transformation 

x — > X — x X G R 

corresponding to the reflection around the point X /2 on the real x axis. This 
kind of transformation will be useful in the study of properties of potentials not 
centred on the origin, and can be readily expressed as the (non commutative) 
product of the parity operator V and the coordinate shift operator V(X ) 

V G (X )=VV(X ). (36) 

Under generalized parity, the kets \R > and \L > are transformed as follows 

V G {X )\R>=e ikX °\L> 
V G {X )\L >= e~ ikX °\R> 

and, correspondingly, the bras < R\ and < L\ give 

< R[P G \X Q ) =< L\e- ikXo 

< LlV^iXo) =< R\e ikx ° 

Since the matrix V G (X ) associated to this transformation reads in the \R >, \L > 
basis 



/ q e -ikX \ 

v g {Xq) = I eikXo I 



and, consequently, 

V G (X ) = V G \X ), 
the four basic operators f^- transformed under V G (X )n i jV G 1 (X ) yield 

&LR &>LL ^RL^~ 2lkX ° QrR 

Thus, the transformed Hamiltonian reads 
H Vg = V g (X )HV g 1 (X ) 

= H RR Q LL + H RL tl LR e 2ikX ° + H LR Q RL e- 2ikX ° + H LL Q RR 

Invariance of the Hamiltonian under generalized parity, H-p G = H, is equivalent 
to 

H RR = H LL (38) 

H RL = H LR e~ 2ikXo . (39) 

The behaviour of the S matrix under V G immediately follows from the V G in- 
variance of H 

S P = V G (X )SV G 1 (X ) 



in 



We have 

S-p G = S , 

or 

V G (X )S = SV G (X ) . 
Vg invariance of the S matrix corresponds to 

Srr = Sn — > Tr->l = T L _> L 

and 

S RL e ikX ° = S LR e- lkXo - R R ^ L e ikXo = R L -> R e 



■ikXo 



3.4 Time Reversal 

For this case we follow the same steps as in subsection 13.21 First, we consider the 
consequence of the invariance property for the Hamiltonian (J26|) . having in mind 
that the time reversal operator T is an antiunitary operator (T = T _1 = T^). 
Under time reversal 

x -> x, p x -> -p x i -> -z 
any ket in the |i2 >, \L> basis is transformed according to 

T(a\R > +(3\L >) = a*\L > +(3*\R >, 

while a bra fulfills the relation 

(7 < R\ + 5 < L\)T~ l = 7* <L\+5* <R\ 

It is worth stressing that T-invariance as such does not force the potential to 
be real, and hence, the Hamiltonian to be hermitian. 
The Qij operators are now transformed into TQy-T -1 

QrL _^ &LL &LR /^q\ 
^LR &LL &RL QrR 

which may be compared to eqs. (|24j) and (|36j) . The T— transformed Hamiltonian, 
Hq- = T7iT~ l , reads 

TT L * L TT* R J ■ ( 41 ) 

RL RR / 

From the properties of the parity operator ([25)1 . it is easy to realize that Hj- 
defined in the previous equation is such that 

VH T = H*P , (42) 



or 



H T = VH*V , (43) 
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where V is the parity operator 

We now proceed to examine how the S matrix transforms under T without im- 
posing any symmetry to the Hamiltonian. Starting from the transition operator, 
T(t — t ), of formula (|55|L and using eq. (jlHj) . we have 

f r (t-t ) = Tfit-t^T- 1 

oo 

= Te iHat T- 1 T{-i) n E n r- 1 (t - t ) n /n\Te' lHata T- 1 

n=0 

oo 

iH t 



e -iH t i n H^(t - t ) n /n\e u 

n=0 

oo 

e -iH t J2 t n (VH*V) n {t - t ) n /n\e lHoto = VT*{t - t )V, 



n=0 

which, when t — > +00 and to ~~ > —00, yields 



S r = TST- 1 = VS*V = ( g L S J; R ) . (44) 

^RL ^RR 



3.4.1 H<r = H 



For a T-invariant Hamiltonian (Hq- = H), the diagonal matrix elements are 
complex conjugate of each other as are the two non-diagonal ones, so that a time 
reversal invariant Hamiltonian is, in that basis, of the form 

H = Ht=( H * R If). (45) 

V H RL H RR J 

Hence, as is well-known, this is not equivalent to hermiticity of the Hamiltonian. 

Coming now to the transition operator, T, we easily obtain, under time re- 
versal invariance of H, with the same procedure as in the previous subsection, 



Hence, H = H? leads to 



T r {t-t ) =T-\t-to) 



1 / Sll —Srl 



Simultaneous validity of Eqs. (}4*4l l4T?j) thus leads to the following explicit 
relations for the S-matrix elements : 

S RL + S* LR detS = 

S LR + S* RL detS = 

Sll = S* LL det S 

Srr = S RR det S . 
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Hence, on the assumption that all the S'-matrix elements are different from zero, 
we obtain 

Sll Srr Srl Slr 



detS 
with 

and 

or, equivalently 



°LL °RR °LR °RL 



| detS |= 1 
Slr 1 = 1 Srl 



I Rl^r \ — \ Rr^l I • 
This also forces SllS rr , i.e. Tr^lT^r to be real and 

L^R-Lr^l+ I J^L^R | — -L , 
R—>L-L L^R+ I J^L^R I — I • 



3.4.2 H T = H ] 

Since, by definition, 



TT* TT* 

tt\ _ I n RR LR 
H RL H LL 



using the explicit representation of Ht, formula (|4Tj) . we may write 

which displays transparently the connection between the time reversal trans- 
formed Hamiltonian and the hermitian conjugate of the same Hamiltonian. They 
coincide if the diagonal matrix elements Hll and Hrr are identical. In that 
case the hermiticity of the Hamiltonian forces time reversal invariance. This is 
what happens for a local potential, yielding the current conservation relation (see 
section HJ. For a non-local potential, in general, Hj- ^ W. 

If we now assume, instead of time reversal invariance {Ht = H), the condition 
Ht = H', equivalent to the intertwining relation 

H T T = TH = H ] T. (48) 

it is almost immediate to check that 

S T = TST- 1 = S f , (49) 



or, in explicit form 



°LL ^LR \ _ / °_RR °LB 



C* c* ; \ Q* C* 
°RL °RR J \ °RL °LL 



leading to the equality of the diagonal matrix elements 

Srr = Sll —> T L ->r, = T R _> L ■ 
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3.4.3 H T = H^ = H 

At this stage, we want to point out that if we assume, in addition to H = Hq-, 
also hermiticity of the Hamiltonian, i.e., H = H>, we see from eq. ()47|) that we 
must have 

TT TT* TT* TT 

n RR — n RR — n LL — U LL ■ 

Thus, the diagonal matrix elements of H are equal and real. 

While condition (J4"H|) may be met by Hamiltonians that are neither hermitian, 
nor T- invariant, for instance those containing local potentials, it is met a fortiori 
by Hamiltonians that are both hermitian and T- invariant. In the latter case, 
we simultaneously have Hq- = H =>- Sq- = S^ 1 , according to Eq. (JSJ), and 
Hq- = W =^ St = S\ on the basis of Eq. (J49j) . From these conditions the 
unitarity of the S matrix follows 

S 1 " = S' 1 . (50) 

In terms of S'-matrix elements, 

S RL + S* LR detS = 

S LR + S* RL detS = 

Sll = S* RR det S 

Srr = S* LL det S 

Summing up, T-invariance and hermiticity of the Hamiltonian lead to 

S l V = VS and S f = S' 1 (51) 



or, explicitly, 



det 5" = — — = -— and S LL = S RR (52) 

^LR ^RL 



so that 



I Slr 1 = 1 Srl I and | det S \— 1 (53) 
In terms of transmission and reflection coefficients, we then have 

Tl-+r = T R -+l , 

I Rl^r I = I Rr-^l I • 

It is worth pointing out that our definition of T is consistent with that of 
ref. while ref. [H| adopts the following 

T = OK, (54) 

where K is the complex conjugation operator, and O is a unitary operator. 
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3.5 VT Symmetry 

We may now try to understand the consequences of PT- invariance (but not 
separately parity or time reversal invariance). 

Any ket in the \R >, \L > basis is transformed under VT according to 

VT(p\R > +a\L >) = p*\R > +a*\L > (55) 

while, correspondingly, for a bra 

(p < R\ +a < L\)T~ X V~ X = p* < R\ + a* < L\ , 

where p and o are complex constants. 

The dyadic operators fly = \i >< j\, = R,L) are thus unchanged under 
VT transformations. Since VT is antilinear, the PT-transformed Hamiltonian 
reads 

H VT = VTHT-^ ^ ^ (5g) 

3.5.1 H VT = H 

According to Eq. (|56|). invariance of the Hamiltonian under VT leads to the 
equality 

Hrr Hr L \ _ f H* RR H* RL 



\ li lr n LL J \ H LR H LL J 

which shows that the matrix elements have to be real. 

From the Schrodinger equation (j3J) and the definition of the L — R basis, it is 
immediate to check that, for a local potential, V(x), 

/+oo 
V{x)dx = H° LL + V LL = H LL . (58) 
-oo 

In order for Hrr = Hll to be real, it is necessary and sufficient that the 
imaginary part of the integral on the r. h. s. of Eq. (J58j) vanishes, i. e. the 
imaginary part of potential V is an odd function of x. As a consequence of the 
reality of the off-diagonal matrix elements of H, i. e. of potential V, one readily 
obtains 



V r (x)sin(2kx)dx = 

/+oo 
Vi(x)cos(2kx)dx = 0. 
-oo 



oo 
oo 



The second relation is automatically fulfilled, since Vi(x) is an odd function of x, 
while the first relation forces V r (x) to be an even function. 
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It is easy to find how the transition operator defined in eq.()3ip is affected by 
VT invariance (reality) of the Hamiltonian, H . Indeed 

oo 

T(t-t ) = Y^(-i) n H n {t-t ) n /n\ 



n=0 

oo 

J2(i) n H n (t-t ) n /nl 

\n=0 

= (T-\t-t )Y, 

and, since 

T- 1 (t - t ) = e^T" 1 (t - t ) e~ lHot , (59) 

T* (t - t ) = e~ iHot T* (t - t ) e iHoto , (60) 
the equality T" 1 (t — t ) = T* (t — t ) implies that 

lim f - l (t - t ) = lim f * (t - t ) , 

t—*-\-oo,to—*— OO t— *+00,t()— > — oo 

hence 

S" 1 = S*. (61) 
This yields for the S -matrix elements 



S^ + S^detS^O 
S L R + S* LR detS = 
S LL = S* RR detS, 
Srr = S* LL det S . 



(62) 



This imposes that S R lS} jR = R r ^lR* l ^ r is real, and that 

| det S \= 1 and | S LL \ = \ S RR | , (63) 
where the latter condition corresponds to 

I T L _ R | = | T R ^ L | . (64) 

3.5.2 H VT = H and H T = W 

Remembering that the time reversed Hamiltonian Hq- is, by definition, such that 

TH T = HT and H T T = TH 

and assuming, in addition, condition certainly valid for a local potential 

and repeated here for clarity's sake, TH = H^T, forces the diagonal elements 



of the Hamiltonian matrix to be equal, Hll — Hrr ; we then see that for a 
PT-invariant Hamiltonian the following relations hold 

VH ] = HP or H*V = VH 

On the conditions given above, 

VHV- 1 = E\ 

i.e., H is pseudo-hermitian [T3j with respect to V. 

As a consequence, the transition operator behaves as follows 

oo 

VT\t-t ) = ^(irP# tn (t-t )7n! 

n=0 
oo 

= ^{i) n H n (t - t ) n V/n\ 

n=0 

= T-\t-t )V. 

Hence, recalling Eq. (J2U), 

VT^(t-t ) = Ve iHoto e iHt (t-to) e -iH t ( 65 j 

_ ■p e iH t 0f J-VHV' 1 {t-t ) e -iH t /gg) 
_ <p2 e iH t 0e iH{t-t ) e -iH tj)-l (g^ 
_ e iH t e iH(t-t ) e -iH t<p (gg^ 

= f- l {t-t Q )V (69) 

and, in the to ~~ > — oo, t — ► +00 limits, 

p^t = 5-1^ (70) 

Thus, under both PT-invariance of the Hamiltonian and the condition that 
H? = H\ so that the diagonal matrix elements of the Hamiltonian are real and 
equal, Hll = Hrr, we have, since S^ 1 = S*, 

VS^ = S*V or S t V = VS. (71) 

In turn, this leads to the equality of the diagonal elements of the S matrix 

Srr = Sll =>- Tl-+r = Tr^l ■ (72) 

Thus 

S RL + S RL detS = 
S LR + S* LR detS = 

Sll = S RR = S* LL det S = S* RR det S , 
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which leads to 

j , o _ Srr _ $ll _ _ Srl _ _ Slr 

c* c* c* c* ' 

J RR °LL °RL °LR 

or, in terms of transmission and reflection coefficients, 
| det^ | = \T L ^ R /T* L ^ R \ = 1, 

dLRORL — dRLJLR - > t^L^RtiR-^L ~ > 

Rl^r + R*l-,rTl^r/TI^ r = , 

Rr^l + R r ^ l Tr->l/T R -^l — • 

It is worthwhile to stress again that the equality (J72*Jl of the two transmission 
coefficients is not a consequence of VT symmetry, which yields only the equality 
of their moduli, but of the additional intertwining condition (|48)L valid for any 
local potential. Again, this point is discussed in Ref. |B]. 

3.5.3 Exact Asymptotic VT Symmetry 



In this sub-section, we consider Hamiltonians with exact VT symmetry, i. e. 
PT-invariant Hamiltonians whose eigenstates are also eigenstates of VT . 

It has been recently proved ^5] that a Hamiltonian with exact VT symme- 
try is unitarily equivalent to a Hamiltonian that is hermitian with respect to a 
suitably defined inner product. 

Let us now investigate some consequences of exact VT symmetry on the S 
matrix. To this aim, it is convenient to introduce the transformation under VT 
of a generic wave function, \I> (x) 

VT^ (x) = ^> V t (x) = (-x) (73) 

and the condition of exact VT symmetry 

^ vr ( x ) = ^* (- x ) = e i0 q> (x) , (74) 

where 6 is a real number, because (VT) 2 = 1. 

Let us apply Eq. (|74^) to the asymptotic wave functions defined in section 2, 

^ vr (±oo) = ^* (too) = e w ty (±oo) , (75) 
and, in terms of their amplitudes, A±, B±, A± and B±, 

A* ± e ikx + B* ± e~ ikx = e ie (A T e ikx + B T e~ ikx ) , (76) 
A* ± e lkx + B* ± e~ ikx = e* (A T e ikx + B T e~ ikx ^j . (77) 

Hence, in particular 

A\ = e w A^ , B* + = e^B. , (78) 
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and 

B_ = e - lQ B\ = , A\ = e ie A_ = . (79) 

The latter equalities come from the boundary conditions for an incident pro- 
gressive wave (B + = 0) and an incident regressive wave (^A^ = J , respectively. 
The S matrix elements are thus written as 

S RR = T L _ R = £ = e~ ie §- = e -*(*+**->, {A _ = e * Q _ } 

SlR = Rl^R = = > 

S LL = T R ^ L = h ; = e-^| = e-<^) , (B + = 
Srl = Rr^l = it = 0. 

The reflection coefficients are thus zero and the transmission coefficients have 
unit modulus, in keeping with the more general condition |det<S| = f, Eq. (|B5|L 
imposed by the VT symmetry of the Hamiltonian. In this case, the S matrix is 
unitary, too, since S^ 1 = S* = (S*) = S^. 

Conversely, it is not difficult to show that conditions (J8Uj) on the S matrix 
elements are sufficient to ensure that the corresponding asymptotic wave functions 
are eigenstates of VT . 

Examples of 'PT-symmetric reflect ionless potentials are discussed in section 

EH 



B 4 



(8Uj 



4 Probability Current and Density for Linear 
and Antilinear Transformations 

In the present section, we consider the time-dependent Schrodinger equation ([TJ 
with a local potential, V(x), and its solution, ip(x,t). In this case, T invariance 
is equivalent to hermiticity and V is real. 

We introduce a linear transformation, Ul, and a corresponding antilinear 
transformation, Ua, commuting with the kinetic energy operator, p 2 = —d 2 /dx 2 , 
and apply them to the Schrodinger equation JTJ with a local potential (jSJ): 

d 2 \ d 

+ v u L (x) i) Uh (x, t) = i^Ul (x, t) ; 



Q X 2 -"I. V~,y -ru L -ft 

d 2 „ , A , , , J) 



- Vu A ( x )J ^u a (x, t) = -i^ipu A ( x i *) 
V v (x) = UV{x)U~ 1 



dx 2 - UA ^') ™A^"J - m 

where = and 



or, equivalently 

V v {x)U = UV{x) , 
namely, if V is invariant under U, it commutes with U. 

1Q 



Now, we multiply the equation satisfied by ipu A by ijj and the initial Schrodinger 
equation by ipu A and subtract them side by side, obtaining: 

(d 2 d 2 \ d 

^ UA dx~^ ~ ^^ Ua ) +{y ~ Vua) M = % di {M) ■ (81) 

Now we introduce the density of probability current: 

ju A (x, t) = -i( (Jl~ij (x, t^j tpu A (x, t)-ij) (x, t) (J^^Ua (x, t) 

and the density of probability: 

Pu A (x, t)=if) (x, t) ip UA (x, t) . 
Therefore, Eq. (j81|) can be rewritten as: 

q^Jua + g- t Pu A =-i(V- V Ua ) PUa . (82) 

If the potential is invariant under Ua the right-hand side of Eq. (|K2|) is zero and 
we obtain a continuity equation, which, for stationary waves, yields a constant 
current. 

Analogously, for a linear transformation, we first have to consider the complex 
conjugate of the equation satisfied by ipu L and repeat the same steps as before, 
by replacing everywhere ipu A by 4>u L and Vjj a by V{} : 

|^ + fp^ = -*(^-^W (83) 

For real V, similar considerations are valid and Eq. (J83j) is reduced to a 
continuity equation. 

The formalism outlined above does not apply to non-local potentials. An 
attempt to extend flux conservation to the latter case is described in Ref . |TB] . 

4.1 The VT- Symmetric Case 

In the case of PT-symmetric potentials, the definitions of p and j satisfying a con- 
tinuity equation are easily obtained by writing the time-dependent Schrodinger 
equations satisfied by ip-p (x, t) and ipr (x, t): 

d 2 \ d 

+ V (-x) ) ip r (x,t) = i—ip P (x,t), 



dx 2 J ' dt 

d 2 A , , .8 



+ V*(x) )ipT(x,t) = -i—^ T (x,t) 



dx 2 J dt 



where ip-p (x, t) = ip (— x, t) and ipr (x, t) = ip* (x, t), so that the equations given 
above are equivalent to the following 

d 2 \ d 

+ V{-x))ip(-x,t) = i—ip(-x,t), (84) 



dx 2 J dt 

<h + v* (x)) r (x, t) = -i^-r (x, t) . (85) 



dx 2 J ' dt 

By multiplying side by side the first equation by ip* (x, t) and the second 
equation by ip (— x,t), and subtracting the two equations side by side, we obtain 

-r (x, t) -f^ip (-x, t) + i> (-x, t) -^r (x, t) + (v {-x) - v* (x)) r ( x , t) ip (-x, t) 

ox dx z 

d d 
= (x, t) —ip {-x, t) + lip (-x, t) —ip* (x, t) . 

VT invariance implies that V (—x) = V* (x), so that the equation above can 
be reduced to a continuity equation 

|W) + |p(M) = o, 

with 

p(x,t)=,l)*{x,t)il>(-x,t) , (86) 

j(x,t) = ±[P(x,t)^(-x,t)-il>(-x,t)^P(x,t)] , (87) 

which is consistent with the definition of j given in Ref. J7j. The following 
relation holds: 

j*(-x,t) = -j (x,t) . 

It is worthwhile to point out that j is identically zero when ip* (x) = e ta ip (—x) , 
where a is a real constant, i. e. ip is an eigenstate of VT. If -0 is a stationary 
wave, ip ( x i t) = e~ lEt ip (x) , j depends only on x and current (|S7jl is constant in 
the whole space. The asymptotic solutions can be taken in the form of stationary 
waves and we obtain, remembering that, for a wave travelling from left to right, 
Eq.©, we have B + = 0: 

If A* + e- ikx (-ik) [A_e~ ikx - B_e +ikx ] 1 
3 [+oo) - - j _ [ A _ e -ikx + B _e +ikx ] (-zife) A* + e~ lkx J 

= 2kA* + B_, 
and, in the same way 

j (-oo) = -2kA + B*_ . 

As expected, j (— oo) = — j* (+oo). Note that, for the reflectionless potentials 
discussed in sub-section 13.5.31 £>_ = . In this case, j (— oo) = j (+oo) = 0. 
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Assuming now B_ ^ 0, conservation of current implies that we also have 
j (—00) = j (+00); hence, the current is purely imaginary and we can write: 



A + B*_ + A* + B_ = , 

or, dividing both sides of the above equation by A_A*_ and remembering the 
definitions of T L ^ R and Rl->r'- 

Rl^rTI^ r + R* L _+ R T L ^ R = , 

a result already obtained in sub-section 13.5.21 or, equivalently: 

Rl->r T l _, r _ 

^L^R 1 L^R 

which means that the phase tp r of the reflection coefficient and <pt of the trans- 
mission coefficient are related by: 

71 

(p r = ip t + — + nir, 

where n is an integer. This relation is trivially checked for square well 
hyperbolic Scarf jH] and generalized Poschl- Teller 0] potentials. 



5 Explicit Relations for Some Potentials 

For the sake of example, the formalism worked out in the previous sections is here 
applied in detail to a few "PT-symmetric solvable potentials. For those already 
considered in Ref. (0); having an imaginary part made of an odd combination 
of Dirac delta functions, QV(x) = A (5 (x + ^f) — 5 (x — ^ L )), it is easy to check 
that all the relations of subsection 13.51 are valid. 

In the presentation of the examples, we make a selection dictated by solvability 
implemented by different methods. 

5.1 Complex VT- Symmetric Square Well 

Let 

V = 0, x < -b , x > +b (88) 

V = -V + iVi , -b < x < , 

V = -V -iV 1 , 0<x<+b, 

Here, V and V are real parameters and V > 0. 

It is worthwhile to mention that discrete states in a PT-symmetric square 
well with Vo = were already studied in Ref.|18j. It was shown, in particular, 
that the square well possesses a real discrete spectrum on condition that the 
coefficient, V\, of the imaginary part is smaller than a certain critical value. On 
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this condition, the eigenfunctions of H are also eigenfunctions of VT, i. e. the 
model has an exact VT symmetry. In this case, it is possible to construct a 
hermitian Hamiltonian unitarily equivalent to the PT-symmetric square well in 
a Hilbert space endowed with a properly defined scalar product, as shown in full 
detail in Ref. [TH]. The analysis has been extended to the continuum of scattering 
states in Ref. [20 J. 
Introducing now 



and 



implying 



with 



k 2 = E 



Or 



Ol 



- E + Vo — iVi 
E + V + iV 1 , 



ao,i = ol exp (Ti<p) , 

Q!q = Oil, a * = a > 



tt' 



{E + Vof + V 2 



iarctanp^ 
2 \E + V 



The two linearly independent solutions (Eqs. (jOllTT|) can both be written in 
the general form 



j^^ikx _|_ j^^—ikx 
Q £ ia x _|_ jj e -ia x 
j^giaxx _|_ p e ~iaix 
Q e ikx + He -ikx 



x < — b 
-b<x<0 
0<x<b 
b < x 



with suitable specification of parameters A, 
The PT-transformed wavefunction reads 



H. 



Q* e lkx + H * e -ikx 
j^* e ia x _|_ p* e ~ia x 
£J* ^ia\x _|_ io.\x 
J^* ^ikx _|_ ~Q* Q—ikx 



X < — b 
-b<x<0 
0<x<b 
b < x 



"PT-symmetric wave functions would thus meet the conditions A = G* , B = 
H* , C = E* and D = F*. In the treatment of scattering states for Vq = 0, 
Ref. |20j considers two linearly independent PT-symmetric wave functions, in 
keeping with the double degeneracy and reality of the eigenvalues E = k 2 . We 
stress, however, that these PT-symmetric eigenfunctions are not the \&i and ^ 
functions (see Eqs. ()12II13|) studied in section IB. 5. 31 



9.3 



We then have three pairs of continuity equations at x — —b, x = and x — b, 
respectively : 




implying 

I C I 2 + I D \ 2 =\ E I 2 + I F 




They allow to re-express the coefficients of the wave function at x = — oo in 
terms of those at x = +00 

/ A \ ( M RR M RL \(G\ 
V B ) \M LR M LL ) \H ) ' 

with 

Mij = M itj (a , a u b, k) = R, L) 

The various matrix elements are given by the following expressions 

M RR = ^e 2ikb {/(«„, on, b, k) + f{a , -ax, b, k) 
+f(-a , ai, b, k) + f(-a , -an, b, k)} 



M L r = g {g(a , an, b, k) + g(a , -a u b, k) 
+g(-a , ai, b, k) + g(-a , -a u b, k)} 



and 



with 



M RL = - {h(a , ai, b, k) + h(a , -a u b, k) 

o 

+h(— ao, ai, b, k) + h(— cto, —a±, b, k)} 
M LL = ^ e ~ 2lkb {m(a , an, b, k) + m(a , -an, b, k) 

o 

+m(—ao, oti, b, k) + m(— a , —on, b, k)} 



f(a ,on,b,k) 



;i + ^)(l + ^)(l + A) e -*(«o+ai)6 
a k oli 



and 

g(a , on, b, k) = (1 + ^)(1 - ^)(1 + A) e -*<°o+«i>* 
/i(e*o, «!, 6, fc) = (1 + -)(1 + ^)(1 - A) e -<(°o+«i)& 

and 

m(ao, ax, b, k) = (1 + ^)(1 - ^)(1 - A) e -*(«o+«i)6 

The obvious symmetries relating these auxiliary functions stem out 

m(a , ai, 6, fc) = f(-a , -ati, -b, k) 

h(a , «i, b, k) = g(-a , -a ± , -b, k) 
From these definitions, we see that 

f(a Q , ai, b, k) = f(ai, a Q , b, k) = f(-a , -ati, -b, -k) = /* a , -b, k) 

h(a ,ai,b,k) = —h(ai,ao,b,—k) = —h(—ai,—a ,—b,k) 

= h(—a , — aii, —b, —k) = —h*(a , oti, —b, —k) 

A series of relations concerning the M-matrix elements follow from the above 
relations: some of which are listed below, not pretending to completeness (note 
that M RR and Mil are invariant under the exchange of ao and a±) : 

M RR (a , an, b, k) = M RR (ai, a , b, k) = M RR (-a , -an, -b, -k) 

= M RR (a ,a!,b,-k) 

and correspondingly for Mll with, in addition, 

M LL (a , oti, b, k) = M RR (a , aii, b, -k) = M RR (-a , -an, -b, k) , 

while 

M RL (a , ai, b, k) = -M RL (ai, a , b, -k) = M RL (-a , -oti, -b, -k) 

= -M RL (a ,a l7 -b,-k) 

and correspondingly for Mlr with, in addition, 

M LR (a , ai, b, k) = M RL (-a , -b, k) = -M RL (ai, a , b, k) 

= -M RL (a ,a l7 -b, k) 

We can then give more explicit expressions for the diagonal matrix elements : 

k 2 - a 2 



M RR = e 2tkb ■ {cos 2 (p cos(2ab cos (p) + sin 2 (p cosh(2ab sin ip) 



2ka 

k 2 + a 2 

sin <p sinh(2a;6 sin <p) — i cos <p sin(2a6 cos ip) 

2ka 



and 

.k 2 -a 2 



2ka 



Mll = e * • |cos y cos(2a&cosy) + sin y cosh(2afesiny) + i- 

k 2 + a 2 

■ sin ip sinh(2a6 sin (p) + i cos (p sin(2afe cos ip) 

2ka 

whereas the non diagonal matrix elements read 

f k 2 - a 2 

Mrl = i { sin <p cos (p[cos(2ab cosy) — cosh(2afe sin ip)] + 



2ka 



k 2 + a 2 

■ cos ip sin(2a6 cos ip) H sin ip> sinh(2a6 sin ip) 

2ka 

and 

f k 2 -a 2 

Mlr = i s sin ip cos (p\cos(2ab cosy) — cosh(2afe sin (p)\ 

[ ' 2ka 

k 2 + a 2 

■ cos tp sin(2a6 cos ip) sin ip sinh(2a6 sin ip) 

2ka 

where we have used the modulus a and the phase <p introduced earlier. 

Transmission and reflection coefficients are easily expressed in terms of the 
M- matrix elements, on the basis of Eq. ( I18|) . 

Assuming an incident wave coming from the left, i.e., from x = — oo, the 
reflection and transmission amplitudes are obtained by assuming H = 

p M LR 



and 



which implies that 



Mrr 
1 



■ L^R — — , 



ill 



RR 



1+ 


M L r 


2 




Mrr 


2 



and 



Rl^r T* l ^r + RUr T l ^r = 2 Re(M LR ) . 
Correspondingly, we have for an incident wave coming from x = +oo 

p _ Mrl 
Rr-^l - -j^ 

and 

_ det M _ 1 

R ^ L ~ Mrr - NIr~r ' 

9fi 



which implies that 



1+ 


M RL 


2 




M RR 


2 



Rr-*l + R* R ^ L T R _+ L = -2 Re(M RL ) . 

The two reflection amplitudes are related by the non-diagonal M-matrix ele- 
ments 

Rr^l = M RL 
Rl^r M lr ' 

It is also easy to check that 

det M = M RR M LL - M RL M LR = 1 , 

which amounts to T L ^ R = T R ^ L . 

It is worthwhile to spend a few words on the behaviour of the transmission 
and reflection coefficients under the exchange of the flux generating part of the 
potential, V(x) = —Vq + iVi, with the flux absorbing part, V(x) = —V$ — %V\. In 
our notations, the exchange V% <-> —V\ is equivalent to ckq «-> a%. Since we have 
already noticed that M RR is symmetric under that exchange, so is Tl->r = 1/M RR 
and, consequently, T r _>l = T^ R . 

The reflection coefficients have a different behaviour: from the relation 

M LR (ai, a , b, k) = -M RL (a , ai, b, k) , 

we immediately deduce that 

Rr-^l(ch , «i, b, k) = R L _+ R (oti, a , b, k) . 

This behaviour of reflection coefficients of PT-symmetric potentials was first 
noticed in Ref. [13J and discussed there in detail for the parallel case of the square 
barrier, as well as of the hyperbolic Scarf potential, to be treated in subsection 

5.2 Multiple Square Well 

In order to extend the formalism of the preceding section to a multiple potential 
well, it is convenient to consider the change of the M matrix elements in the case 
the single well is not centred on the origin, but on an arbitrary point X on the 
real axis. Under the coordinate shift x' = x — X , the basic vectors are changed 
as follows 

= e -ikX |^ . = e ikX ^ ^ 
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and formula (|23|) holds, with a shift C = — X , so that the diagonal elements of 
the M matrix remain unchanged, while the off-diagonal ones are changed by a 
phase factor 

M' R r M' RL \ = f Mrr e~ 2ikX °M RL 
Kr Kl J V e 2tkX °M LR M LL 

We are thus ready to solve the simplest problem of multiple wells, i. e. two 
identical square wells of width 2b separated by a zero-potential region of width 2a 
centered on the origin. The M matrix of this problem is nothing but the product 
of the M matrices of the two wells, with off-diagonal elements properly changed 
in phase so as to take into account the shifts of the well centres with respect to 
the origin 

Mrr e 2tk ^M RL \ ( Mrr e^^Mjn 

e -2ik(a+b) MLR Mn J V e +2 ^ a+b ^M LR Mn 

(89) 

The Mij elements are, of course, the same as in the preceding section. 

In order to extend the above formalism to an arbitrary number, n, of iden- 
tical square wells of width 2b separated by intervals of constant length 2a, it is 
convenient to introduce a new transfer matrix, T, connected with M as follows 



M (dw) = M W. M (2) 



Trr Trl \ = f M RR e- 2lk ^ M RL e 2lka 
Tir T LL ) V M LR e- 2ika M LL e 2ik ^ 



(90) 



It is now easy to check that the M^ matrices of formula (J89*j) can be written 
in terms of T in the following way 

M (i) = V *(-a -2b) -T- V{a), M (2) = V*(a) -T-V(a + 2(a + b)) , 

where T>(x) is the diagonal matrix defined by formula It is worth pointing 
out that the argument u\ = —a — 2b of the first T>* matrix is the coordinate of the 
left edge of the first well, while the argument v i = U\ + 2 (a + b) = a of the first V 
matrix corresponds to the left edge of the second well, obtained by summing to 
the left edge of the first well the spatial period 2(a + b). A similar interpretation 
holds for the arguments Ui = a and vi = a + 2 (a + b) of the second well. 

Therefore, the M matrix of the double square well can be written in the very 
simple form 

M = V*(-a-2b)-T-V(a)-V*(a)-T-V(a+2(a+b)) = V*(-a-2b)-T 2 -V(a+2(a+b)) . 

The generalization to an arbitrary number, n , of square wells with space period 
s = 2(a + b) over an interval of length L = ns is trivial 

M = V*( Ul ) -T n -Viux+ns) . 

In the n — ► oo limit, the above formula corresponds to the PT-symmetric 
version of a model of one-dimensional crystal, and might be used in an analysis 
similar to that already carried out in Ref. [21] in the hermitian case. 



5.3 The Hyperbolic Scarf Potential 



As an another example of solvable potential, we consider the hyperbolic Scarf (or 
Scarf II) potential, whose scattering solutions were investigated in Ref. m the 
hermitian case and in Ref. [3] in the PT-symmetric case. This potential allows 
simple analytic solutions for the transmission and reflection coefficients, thus 
displaying explicitly the singularity structure in the complex k plane characteristic 
of a PT-symmetric potential. General comments on these singularities as poles 
of the S matrix connected with bound states and resonances were recently made 
in ref. We repeat here the hermitian case of the Scarf II potential, because 
the solutions given in Ref. J221 contain several misprints. Following the notations 
of Ref. |H] , the hermitian Scarf II potential is written in the form 

V (x) = (A 2 - s (s + 1)) — L- + X(2s + 1) , (91) 

cosh x cosh x 

where A and s are real parameters. Two independent scattering solutions are: 



F 1 (x) = (l + iy) 2 (1-iy) 2 F -s - ik, -s + ik,i\ - s + -; 




s + l-iX s + iX ( 1 1 3 1 + III 

F 2 (x) = (l + iy) 2 {l-iy) 2 F iX-ik, iX + ik,s-\ zA; 



2 7 2 2 ' 2 

(93) 

where y = sinh x and F (a, b, c; t) is the hypergeometric function. By exploiting 
the following asymptotic formula of the hypergeometric function 23 

, Hm F (a, i, „ t) = r (c) ( r .^7 a) (-tr + ^\ (-*)"' 
|t|-+oo \T(b)T(c — a) r (a) F (c — b) 

and the elementary limit lim^^-i-oo sinh x = j- fE£gg) ; we readily obtain the be- 
haviour of Fi and F 2 at x — ■> ±oo : 



lim Fx (x) = a 1+ e ikx + b 1+ e ikx , 



with 



(95) 



with 



e -f {X-k+is) 


r (ix - s + i) r (2zife) 


2-s+2ik 


r (-s + ik) r (| + iA + ik) 


e -f (A+fc+is) 


r (iA - s + 1) r (-2zife) 




r(-s -iife)r (| + iA -ifc) 


lim Fi (a;) 


= ax-e ikx + b X -e~ ikx , 


gf (A+A:+is) 


r [ix - s + 1) r (-2ik) 



ai ~ 2*-™ T(-s-ik)T(\ + iX-iky (96) 



2Q 



Therefore, the following relations hold 

a = e n(X+k+is) b 



e n{x - k+is) a 1+ . 



The procedure is repeated for the second solution, F 2 , with the following 
results: 



lim F 2 (x) = a 2+ e lkx + b 2+ e 

x— >+oo 



-ikx 



where 



e 2 



(A+fc+i(s+l)) 



o 2 h 



r (s + § - iA) r (2ik) 



^2+ 



2 2ife+iA-i r Q - iA + iife) r (s + 1 + life) 

el (A-fc+*(«+i)) r (s + § - iA) r (-2ik) 



2-'"'-'^-a r (i 



-2ifc+iA- 

lim F 2 (x) 



2 

a 2 _e 



iA 



zfc) r ( s + 1 

2-e" 



ifc) 



(98) 
(99) 



kx 1 ^ „—ikx 



where 



e -§(A-fc+i(s+l)) 



a 2 - 



r (s + § - zA) r (-2iife) 



-2ik+i\- 



e 2 



(A+A+i(s+l)) 



r (| - iA - life) r (s + 1 - life) 

r (s + § - iA) r (2ik) 



Therefore: 



)2ifc+iA-7 



12- 

b 2 - 



T (± - i\ + ik)T (s + 1 + ik) 



e -^x-k+i(s+i)) b2+ . 
e _, (A+fc+i(s+1) ) a2+ _ 



(100) 
(101) 



It is worthwhile to stress that the relations connecting aj_ with b i+ and 6j_ 
with a i+ (i = 1,2) are valid not only for the hermitian potential (real s and 
A), but also for the "PT-symmetric potential with real s and imaginary A = iA', 
provided the x coordinate is real. 

After some manipulations of the V functions in the asymptotic amplitudes, 
it is not difficult to obtain the compact form of the T L ^ R and Rl^r coefficients 
first derived in Ref.f22j, corrected in Ref.jH] for a wrong sign, and repeated here 
for the sake of completeness: 



Rl^r 



T (s - ik) T (s + 1 - ik) r (| + iA - ik) V (\ - iA - ik) 



T, 



L-*R 



r (—ik) r (1 — ik) (r(±-ik)) 

cos (its) sinh (ttX) .sin (its) cosh (ir\) 



cosh (irk) 



sinh (irk) 



; (102) 

(103) 



an 



Formulae ([10 2)) hold also for the "PT-symmetric version of the potential with 
real s and imaginary A = iX' . In this case the hyperbolic functions of A are 
changed into circular functions of A': sinh(iA') = zsin(A'), cosh(zA') = cos(A'). 
Formulae (J102j) satisfy the unitarity condition 

\T L ^ R \ 2 + \R L ^ R \ 2 = 1 . (104) 

in the hermitian case (real A), while in the PT-symmetric case (imaginary A) 
unitarity may be broken : for instance, when s is integer and X/i half- integer, 
\Tl^r\ 2 + |-Rl^_r| 2 —>■ oo when k — > . 

By using the general definitions of transmission and reflection coefficients, 
considered as functions of the coupling strength A, together with the relations 
connecting aj± and bi± {i = 1,2), it is easy to check that Tr_^l(A) = Tl^r(X) 
and Rr-^l(X) = Rl^r(-X) in both the hermitian and the PT-symmetric case 
mentioned above, as already noticed in Ref . [T3] . 

A further transformation preserving PT symmetry of the Scarf potential with 
real s and imaginary A is the complex coordinate shift x — > x + ie (— | < e < +f , 
in order to avoid singularities in the potential). 

The asymptotic forms of the two independent solutions F\ and F 2 are easily 
computed by the procedure described in the previous pages: 



s — iX . k s-\-iX 

lim Fi (x + ie) = lim I i I ( —i- 



x^+oo x— >+oo 



F I — s — ik, —s + ik,iX — s + -; i 



2' 4 

a 1+ (e) e ikx + b 1+ (e) e~ ikx 



where 



ai+ (e) = a 1+ e~ ke ; (105) 



b 1+ {e) = b l+ e k \ (106) 

Here, a i+ and 6 1+ on the right-hand-side of the previous equations are given 
by formulae and ([§5]). respectively. In the same way, one gets 



lim F 1 (x + ie) = a x _ (e) e lKx + 6i_ (e) e~ lkx , 

x^—oo 

with 

ai_ (e) = a^e~ ke ; (107) 

b x _ (e) = b^e ke . (108) 

Here, ai_ and &i_ on the r.h.s. are obviously given by formulae 1)96)1 and (|97j). 
respectively. 
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The two limits of F 2 (x + ie) are computed in the same way 



lim F 2 {x + ie) = a 2± (e) e ikx + b 2 ± (e) e~ ikx 

a;— >±oo 



with similar results: 

a 2+ (e) = a 2+ e~ kt ; (109) 

b 2+ (e) = b 2+ e ke - (110) 

a 2 _ (e) = a 2 _e" fce ; (111) 

&2-00 = & 2 -e fc£ . (112) 
It is thus immediate to check the result of Ref. [3]: 

T L ^ R (e,X)=T L ^ R (0,X); (113) 



Rl^r (e, A) = R L ^ R (0, X)e 2ke , (114) 

where the Tl^ r and Rl-*r coefficients on the r.h.s. of formulae (jll3H114|) are 
given by (|102|) . Unitarity is obviously broken by e ^ 0. Moreover, it easy to 
check that, in the same case 

T R _> L (e,X) = T L _> R (e,X) ; (115) 



Rr—*l (e, A) = R L ^ R (e, -A) e" 4fce . (116) 
5.4 Reflect ionless Potentials 

Hermitian reflectionless potentials are much studied in the literature |24[ 1251 12b] . 
In this section, we discuss two examples of reflectionless PT-symmetric potentials. 
The first example is the regularized one- dimensional form of the "centrifugal" 
potential 

V(x) = -^— 2 , (117) 

[X + IE) 

where a is a real strength and e a real constant that removes the singularity at 
the origin. The time-independent Schrodinger equation for the potential under 
investigation reads, in units % = 2m = 1 

d 2 a \ o , 

+ 7 "I J * = k ® > (H8) 



dx 2 ( x + ie)" 
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We introduce the complex variable z = k (x + ie) and express Eq. (|118|) in 
terms of z, 

n d 2 



dz 2 



V + (z 2 - a) * = 



;ii9) 



If we define the new function $ (z) such that \I/ (#) = z 1 ^ 2 ^ (z), the equation 
fulfilled by $ is promptly obtained from Eq. ()119p in the form 



cP 



d 



Z --t;® + Z — $ + Z 



dz 2 



dz 



a: 



(120) 



i.e. a Bessel equation with square index v 2 = a + 1/4. 

A couple of linearly independent solutions to Eq. (jl20)l with the appropriate 
asymptotic behaviour for ^ to be a scattering solution to Eq. (jll9|) is provided by 
the Hankel functions of first and second type, (z) and (z), respectively, 
whose lowest order asymptotic expansions are |2Zj 



lim H<P(z) 
lim Hl 2) {z) 

\z\— >oo 



1/2 

— J exp 

1TZ 



TV 
— l 

2 



(— 

\1TZ 



1/2 



exp 



7T 
— i 

2 



t)1 




4/J 




7T 




~ 4 


)] 



121) 

;i22) 



valid for 3ftz/ > —1/2, |argz| < n. 

The corresponding asymptotic solutions to Eq. ([119)1 thus are 

lim (x) = exp ( ikx — ke — i—v — i— ) , 
x-^oo V 2 4 / 

/ 7T 7T \ 

lim \I/ 2 (a?) = exp I —ikx + ke + z— i/ + z— - I 

x^oo V 2 4/ 

If the above asymptotic wave functions are written in the general form 



lim (x) = eij± exp (z'fcx) + &y- exp (—ikx) , 

x— >±oo 



we immediately obtain 



a 1+ 

«2+ 



= exp 
a 2 _ = ; 6 2 



fee — z— ^ 
2 



7T 



h+ = h- = 



6 2 - = exp ( +fe + z-z/ , , ( 



(123) 
(124) 

(125) 

(126) 
(127) 



Formulae (|12fij) show that VT is an exact symmetry for this potential. In fact, 
remembering Eqs. (j78H79jl . we have in this case A + = ai + , A- = a\- = A + , so 
that A* + /A_ = a\ + /ai + = exp (in (v + 1/2)). In the same way, we get B + = b 2 +, 
B_ = b 2 _ = B + , and B% /B_ = b* 2+ /b 2+ = exp (-in (v + 1/2)). 
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The resulting transmission and reflection coefficients for waves travelling from 
left to right and vice versa are promptly evaluated from their definition 



T^ R = "^ =1, (128) 

a 2 -Oi + — ai-6 2+ 

= 61+62- - 6!_6 2+ = 
a 2 -Oi + — ai-0 2+ 

= g^u. , (130) 

a 2 _6i+ - ai_6 2+ 

= ai+ T - ai t 2+ = • ( 131 ) 

The second example is the hyperbolic Scarf potential with integer coupling 
strengths, s = n and A = im . The formulae of the preceding section read in this 
case 

T - t _ / 1 ,n+m (" ~ -ik)(m-\- ik) ... (| - ifc) ^ 

L ^ " ^ L_( j (n + ^)...(l + zfc)(m-| + zA;)...(| + zA;) (!ld2j 
i^_> fl = i2 M = 0. (133) 

with IT^J = 1. Equations (|132H133j) are sufficient to ensure that the asymptotic 
wave functions are eigenstates of VT. 

Moreover, the potential possesses bound states, which, in the present parametriza- 
tion, turn out to be eigenf unctions of VT , as discussed in Refs. [23 I2H] • 

In the Hermitian case, the hyperbolic Scarf potential is reflectionless only 
when s = n and A = 0, corresponding to the well-known case of the Poschl- Teller 
potential with integer coupling strength. 

5.5 A Non-Local Potential 

Let us go back to the general Schrodinger equation (JTJ for a monochromatic wave 
(J2J) of energy E = k 2 



d 2 



dx 2 



V(x) + A j K(x, y)^(y)dy = k 2 ^(x) , (134) 



where the potential strength, A, is a real number. It is easy to check, by calcu- 
lating scalar products, that the kernel of a hermitian non-local potential satisfies 
the condition 

K(x,y) = K*(y,x). (135) 

In the L — R basis, K is written 2 hermitian matrix, as a consequence 

of the above constraint ()135j) 

K = ( Krr Krl \ = ( K * RR K * LR \ = K-t 
- Kut K LL ) \ K* RL K* LL ' 
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Parity invariance of the potential could be similarly checked to imply 

K(x,y)=K(-x,-y) . (136) 
In the L — R basis, this corresponds to 

r* L )=( f LL t LR 1 = VKV~ l = K v , (137) 

V K LR J^LL J \ J^RL J^RR J 

as in Section l3~2l 

The condition of time reversal invariance of K is obtained by using the defi- 
nitions of Section 13.4.11 in the form 



K T = TKT~ l = VK*V =K , (138) 



or, in matrix notation, 



Krr K rl \_ K* ll K* lr 



K LR Kll J \ K* RL K* RR 



(139) 



which corresponds to 

K(x,y) = K*(x,y). (140) 

Conditions similar to those of Section 1331 therefore, lead us to introduce VT 
invariance of K in the form 

K„ S VTKT-^ = ( K kTr K K f L ) = ( f - K £ l )=K. (14!) 
This corresponds to 

K(x,y)=K*(-x,-y) , (142) 

in agreement with formula (3) of Ref. jHj , which corrects a misprint in the corre- 
sponding formula (113) of Ref. jH]. 

In order to deal with a solvable PT-symmetric potential, we consider only 
separable kernels of the kind 

K(x,y)=g(x)e ia *h(y)e i Py, (143) 

where a and (3 are real numbers, and g(x) and h(y) are real functions of their 
argument, suitably vanishing at ±oo. 

For this kind of kernel, the hermiticity condition ()135|) implies a = —f3 and 
g = h. Parity invariance ()136|) requires a = (3 = and g(x) = g (—x), h (x) = 
h (—x). Time reversal invariance (|14Uj) requires a = (3 = , but does not impose 
conditions on g and h. 

The various conditions that can be imposed on kernel ()143|) are summarized 
in Table I. Finally, VT invariance ()142j) does not impose conditions on a and 
j3, but requires g(x) = g(—x),h(y) = h(—y). As an important consequence, 

3F> 



Reality 


a = (3 = 


Symmetry under x <-> y 


a = (3, g = h 


Hermiticity 


a = —(3, g = h 


V Invariance 


oc = (3 = 0, g (x) = g (-x), h(y) = h (-y) 


T Invariance 


a = (3 = 


VT Invariance 


g(x) =g(-x), h(y) = h(-y) 



Table 1: Possible symmetries of the separable kernel ()143j) 

their Fourier transforms, g(q) and h(q'), are real even functions, too. In order 
to solve eq. (|134j) . we resort to the Green's function method. As is known, the 
Green's function of the problem is a solution to Eq. (|134|) with the potential term 
replaced with a Dirac delta function 

-^G ± (x, y) + (k 2 ± ie)G ± (x, y) = 5{x - y) . (144) 

Here, we introduce the infinitesimal positive number e in order to shift up- 
wards, or downwards in the complex momentum plane the singularities of the 
Fourier transform of the Green's function, G±(q, q'), lying on the real axis. 

In fact, after defining the Fourier transform, f(q), of a generic function f(x) 

as 

/+oo 1 p+oo 

f(x)e-^dx <-> /(*) = —/ f(q)e^dq, 

and expressing G±(x,y) and 8{x — y) in terms of their Fourier transforms, we 
quickly solve eq. ()144|) for G± 

~ 2n5(q + q') 

G±{q,q) - 



—q 2 + k 2 ±ie 

Therefore, the Green's function in coordinate space is 

G±{x,y) = -— / 2 r . e^-^dq = G ± (x - y) . (145) 
2tt J^oo q l - k l =F is 

The integral ()145|) is easily computed by the method of residues. In fact, the 
integrand in G + (x — y) has two first order poles at q\ = k + ie' and q 2 = —k —ie', 
where e' = e/(2k): the integral is thus computed along a contour made of the real 
q axis and of a half-circle of infinite radius in the upper half-plane for x — y > 0, 
on which the integrand vanishes, thus enclosing the pole at q = qi, and in the 
lower half-plane for x — y < 0, enclosing the pole at q = q 2 , for the same reason. 
Notice that the G + contour integral is done in the counterclockwise direction for 
x — y > 0, while it is done in the clockwise direction for x — y < 0, so that the 
latter acquires a global sign opposite to the former. The result is 

G + {x -y) = ~ [e ik{x - y) 6{x - y) + e~ ik{x - y) e{y - x)] , (146) 
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where 9(x) is the step function, equal to 1 for x > and otherwise. 

The second Green's function, G^(x—y), is the complex conjugate of G + (x — y) 

G_(x-y) = ^- [e- lk{x - y) 6{x -y) + e ik{x - y) 6(y - x)] . (147) 

Now, we go back to eq. (j!34|) with kernel (|143|) . call *&±(x) two linearly 
independent solutions, for a reason that will become clear in the next few lines, 
and define the following integral depending on ty± 



-oo 



I ± (J3,k)= / e^h(y)^ ± (y)dy. 



It is easy to show that I±((3,k) can be written as a convolution of the Fourier 
transforms of h(y) and ty±(y). 

The general solution to eq. ()134|) is thus implicitly written as 

/+oo 
G ± (x-y)g(y)e iay dy. (148) 
■oo 

Eq. ()148|) allows us to express I±(/3, k) in terms of the constants c± and d± 
and of Fourier transforms of known functions : in fact, by multiplying both sides 
by h(x)e lf3x and integrating over x, we obtain 

I±(J3, k) = c±h(k + P) + d±h(k -(3) + \N±(a } (3, k)I±(j3, k) , (149) 

where we have exploited the symmetry h(—k — (3) = h(k + (3) and N± is defined 

as 



N±{a,f3,k) = / h{x)e il3x G ± {x -y)g{y)e iav dxdy (150) 

J — oo 

= T —j h{x)e ipx e ±ik \ x - y \g{y)e iay dxdy, 



so that 



= c ± h(k + p) + d ± h{k (3) = ~ + d± ~ h{k _ 
1 — Al\±(a, p, k) 

where 

D ± (a,/3,k)= ' 



1 - \N±(a,/3, k) ' 

Let us examine now the asymptotic behaviour of the two independent solu- 
tions, starting from ^+{x) 

/+oo 
G + (x-y)g{y)e iay dy. (152) 
■oo 
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The asymptotic behaviour of the integral on the r. h. s. of eq. (|152|) is promptly 
evaluated by observing that, according to eq. (|146|) . 



lim G + (x-y) 

x—>±oo 



e ±ik(x-y) 

2k 



so that 



lim \Jf+(x) = c + e lkx + d+e %kx - iul+(ft, k)g(k =F a)e 

x— >±oo 



±ikx 



where we have put u = \/{2k). 

Remembering the expression (jl51j) of J + , we finally obtain 



lim *f?+(x) = c + e ikx + \d + - iug(k + a) c + h(k + ft) + d+h(k - ft) 



Da 



lim ^A r (x) 

x— >+oo 



c + — iujg{k — a) c + h(k + ft) + d + h(k — ft) 



Da 



Akx 



d+e 



-ikx 



ikx 



The constants c + and d + are fixed by initial conditions: if we impose that 
^+{x) represents a wave travelling from left to right, according to formula (|12p. 
we immediately have c+ = 1, d + = and 



Tl-^r = 1 — iujg(k — a)h(k + ft)D + (a, ft, k) 
Rl^r. = —iug{k + a)h(k + ft)D + (a, ft,k) . 



(153) 



It is worthwhile to point out that the above expressions break unitarity, i.e. 
T L ^ R \ 2 + | Rl^r | 2 t^ 1, because probability flux is not conserved in general. 
We come now to the second solution, \&_(a;), written in the form 



+oo 



*_(a:) = c_e lkx + d.e~ lkx + XI. (ft, k) / G_(ar - y)g(y)e iay dy . (154) 



The asymptotic behaviour of the Green's function, G-(x), is now 



Azik(x-y) 



so that 



lim G_(x — y) = — e 
x^±oo 2k 



lim = c-e lkx + d.e~ ikx + iwI-(P, jfe)^(ife ± a)e Tifcc 

x— >±oo 



or, using the explicit expression (j!51j) of J_, 

lim = d_e _ifcB + jc_ + - a) [c_/i(Jfe + /3) + dSh(k - ft) 



lim = c„e ite + <! d_ + iwc7(Jfe + a) c.h(k + ft) + d_h(k - ft) 

x— >+oo 



MS 



Since and *&+(x) are linearly independent, we can impose that 

is a wave travelling from right to left, according to formula the initial con- 
ditions are 

c^ + iug(k-a)(c-h(k + (3) + d-h(k- P))D_(a,P,k) = 0, 
d_ + iLug(k + a)(c_h(k + f3) + d„h(k - f3)))D_(a, f3,k)) = 1, 

where 

d- = T R ^ L , c_ = Rr-, l • (155) 

We then obtain 

Tr^ l = 1 - iug(k + a)h{k - f3)V_(a, p,k) , (156) 
Rr-^l = —iujg(k — a)h(k — /3)V_(a, j3, k) . 
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where 

T>(a,(3,k) = 

1 - AiV_ + za;(<7(& + a)/i(Jfe - (3) + ~g{k - a)h(k + /3)) 

Formulae (jl54H155|) show that, in general, for a PT-symmetric non-local po- 
tential, Tr_,l 7^ T L ^ R . In fact, from the quoted formulae, 

Tr^ l - T L ^ R = iuAD + (a, (3, fc)ZL(a, fc) , 

where 

A = g(k - a)h(k + (3) - g(k + a)h(k - f3) 

+X(N + g(k + a)7i(A; - /3) - N-g(k - a)h(k + (3)) 

+iug(k - a)h(k + (3)(g(k + a)h(k - (3) + g(k - a)h(k + (3)) . 

Computation of the N± integrals yields the general forms 
XN+(a,f3,k) = -%- g(k - a)h(k + (3) + g{k + a)h(k - (3) +Q(a,(3,k), 

XN- (a, /3, k) = i^ g(k - a)h{k + (3) + ~g(k + a)h(k - (3) +Q(a,(3,k), 

where the function Q (a, f3, k) is real. 

If we now make the additional assumption that our kernel is symmetric, 
K(x,y) = K(y,x), i.e. g = h and a = (3, we obtain T R ^ L = T L ^ R . It is 
worthwhile to stress that imposing the symmetry of the kernel is equivalent to 
imposing the intertwining condition (jlSj) . i.e. Kq- = K\ which ensures the equal- 
ity of the two transmission coefficients. 

A detailed calculation for the case 

g(x) = e-rW , %) = e-*l"l , 

with 7 and 5 positive numbers, i.e. the Yamaguchi potential, has been presented 
in Ref. [U]. 

Unitarity properties are discussed therein: in particular, when the Yamaguchi 
potential is PT-symmetric with non-zero a and f3, the characteristics of unitarity 
breaking are distinctly different from those of PT-symmetric local potentials 
discussed in Ref. [T3] . 
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6 Conclusions 



In this final section, we try to focus attention on what we believe are the most 
original results of our investigation. 

Exact VT invariance has been introduced in the literature as a condition on 
the bound-state eigenfunctions of a VT -symmetric Hamiltonian. The conclusion 
is the following: the eigenstates of H should be eigenstates of VT , too. In turn, 
this condition ensures that the corresponding eigenvalues are real. 

As a particular example, we mention the imaginary PT-symmetric square well 
studied in Ref . [TH] . There, a well defined threshold was found for the discrete spec- 
trum, separating the regime of exact VT symmetry from that of spontaneously 
broken symmetry. In the case of scattering, an extension of these results might be 
ambiguous, in so far as the continuum can always be labelled with a real energy , 
E = k 2 . Correspondingly, Ref.[2U] shows that one can always find PT-symmetric 
continuum eigenfunctions of this Hamiltonian, both below and above the critical 
potential strength found in Ref. JXHI • 

In our considerations, we have tried to introduce a specific "exact" VT invari- 
ance associated with the scattering states of type (|T2|) and (|T3J). We have called 
this this condition exact asymptotic VT symmetry, which is only and specifically 
relevant to scattering; we have shown that this condition forces the PT-symmetric 
potential to be reflectionless ( the above mentioned PT-symmetric square well 
does not belong to this class). 

While the interest in reflectionless potentials was recently revived in the frame 
of supersymmetric quantum mechanics and Darboux transformations [26J, we 
stress the fact that the potentials considered in Ref.|2Sl are real. Our link of 
exact asymptotic VT invariance with reflectionless complex PT-symmetric po- 
tentials should provide the necessary stimulus to broaden the investigational ITU] 
and classification of reflectionless potentials so as to include their complex form, 
of which we have provided a few examples in the present work. 

The other most important topic we have elaborated concerns the delicate 
distinction and interplay between hermiticity and time reversal invariance for 
non-local potentials and a proper extension of VT invariance to this case. 

An explicit construction of a solvable separable complex potential has been 
presented and worked out in detail in Ref. [3U]. A particularly notable difference 
between local and non-local PT-symmetric potentials is the non-equality of the 
two transmission coefficients in the non-local case; they have also a quite different 
behaviour in unitarity breaking [P31 IHO] . 

An extension of the present work can be envisaged for multi-channel problems [311 
021 EB] , with the specific aim to formulate a self-contained and consistent frame- 
work to extend the discussion of symmetry properties to elastic scattering of 
non-zero spin particles and to inelastic scattering. 
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